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ABSTRACT 

To assess the economic impact of bovine coronavirus 
(BCV) and bovine respiratory syncytial vims (BRSV) 
infections, accurate estimates of their associated effects 
on animal performance are needed. This study aimed to 
quantify the variation in individual test-day milk yield 
and somatic cell count, risk of reproductive failure after 
first service of dairy cows, and risk of death of calves 
and heifers according to the BCV and BRSV status of 
the herd. Three types of status were defined for BCV 
and BRSV infections, based on 1) the dynamics over a 
7-mo period of BCV- and BRSV-specific antibody levels 
in pooled milk of primiparons cows; 2) the possible oc¬ 
currence of presumably BCV- and BRSV-related clini¬ 
cal outbreaks; and 3) the combination of both pieces 
of information. A total of 36,184 test days, 2,716 cows 
with a first service, and 4,104 calves and heifers in 65 
Swedish herds were included in the analyses. Animal 
performance associated with BCV and BRSV infections 
was quantified using hierarchical mixed generalized and 
survival models, after adjustment for covariates known 
to influence the performance under study. A significant 
reduction in milk yield was observed for cows in presum¬ 
ably BRSV recently infected herds, as well as in herds 
having a presumably BRSV-related clinical outbreak 
(of 0.57 and 0.91 kg/d, respectively), compared with 
cows in presumably infect ion-free herds. There was also 
a significant increase in somatic cell count (of 12,000 
cells/mL) for cows located in herds with a BRSV out¬ 
break. The risk of failure after first service, as well as 
the risk of death in calf and heifer, was numerically 
higher in BRSV-infected herds, although this was not 
statistically significant. In contrast, BCV infection 
herd status, as defined in the present study, was not 
significantly associated with any production losses in 
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animals from infected herds compared with those in 
infection-free herds. 

Key words: dairy cattle, animal performance, bovine 
coronavirus, bovine respiratory syncytial virus 

INTRODUCTION 

Bovine coronavirus (BCV) and bovine respiratory 
syncytial virus (BRSV) are two viruses with worldwide 
distribution that are frequently involved in the respira¬ 
tory and enteric disease complexes of cattle (Clark, 
1993; Valarcher and Taylor, 2007). Bovine coronavirus 
has been described to be involved in winter dysentery 
in adults (Saif, 1990; Alenins et ah, 1991), calf diarrhea 
(Stair et ah, 1972), and in respiratory disease of young 
stock (McNulty et ah, 1984). Bovine respiratory syn¬ 
cytial virus is recognized as one of the most important 
causes of respiratory tract disease in beef and dairy 
cattle, especially in young animals (Stott et ah, 1980; 
Verhoeff et ah, 1984; Elvander, 1996). 

These infections are widespread in dairy cattle. 
Among 113 dairy herds in the Ithaca area (New York, 
NY), a 33% incidence of winter dysentery was recorded 
over a 9-mo period (White et ah, 1989). Bulk tank milk 
data, reflecting the antibody status of lactating cows, 
also showed very high nationwide prevalences of BCV 
(89.2 to 100%; Paton et ah, 1998; Traven et ah, 1999) 
and BRSV infections (70 to 100%; Elvander, 1996; Pa¬ 
ton et ah, 1998) in the UK and Sweden. The use of a 
BRSV vaccine is considered to have a minimal effect 
on antibody levels, as its implementation was rare in 
UK (Paton et ah, 1998) and not used at all in Sweden. 
Furthermore, there is no BCV vaccination program in 
Sweden, and none was mentioned in Paton et al. (1998). 
Moreover, a high BCV and BRSV infection burden and 
spread in calves has been recently demonstrated in the 
southwestern part of Sweden (Hagglund et ah, 2006). 
Measures to control the spread of these infections may 
be warranted, because infectious diseases are likely 
to have detrimental effects on the animals and their 
productivity. However, information on the production 
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effects associated with these infections on a dairy farm 
remains scarce, and an assessment of production effects 
at the farm level might contribute to support decision¬ 
making for control. 

Practical experiences of farmers indicate that the 
losses associated with these infections may be substan¬ 
tial because of reduced productivity. Increased risk of 
mortality and premature culling, reduced growth, fer¬ 
tility, and milk production in first lactation associated 
with bovine respiratory and enteric disease in dairy 
calves and heifers have been quantified in the literature 
(Waltner-Toews et ah, 1986; Curtis et ah, 1988, 1989; 
Warnick et ah, 1994, 1997; Virtala et ah, 1996). How¬ 
ever, these studies did not focus specifically on BCV or 
BRSV infections or on adult cows, which constitute the 
bulk of the production animals in a dairy herd. 

For BCV infection specifically, mortality of calves in 
some case herds with acute outbreaks has been described 
(Reynolds et ah, 1986), and a decrease in milk yield in 
relation to an extended decrease in appetite was also 
reported in both experimental (Traven et ah, 2001) and 
field (Saif et ah, 1988) studies about winter dysentery 
(associated with BCV) in lactating cows. Some severe 
decreases in milk yield associated with a BRSV out¬ 
break have been described by Odegaard and Krogsrnd 
(1977) and Harrison and Pnrsell (1985). Only one study 
(Norstrom et ah, 2001) has quantified this loss, with ac¬ 
counting for other factors affecting milk yield: the daily 
loss was estimated to be small (0.7 kg per cow and day 
in the week following a BRSV outbreak, compared with 
the period >1 wk before the outbreak), but the authors 
stated that it was probably underestimated. 

To our knowledge, no previous investigations have 
focused on the possible variation in performance of 
cows and young stock in commercial dairy herds with 
different statuses of BRSV and BCV infections. There¬ 
fore, the aim of the present study was to quantify, in a 
sample of Swedish dairy farms, the possible variation in 
milk yield, SCC, and reproductive performance of cows, 
and in the risk of death of young stock, according to the 
status toward BCV and BRSV infections of the herds. 
Our assumption was that a recent or current BCV or 
BRSV infection could be associated with a reduction 
in individual milk yield, higher SCC, and reduction in 
reproductive performance (in relation to fever and loss 
of appetite) of cows and a higher risk of death of young 
stock. 

MATERIALS AND METHODS 
Herds and Animals 

Data came from dairy herds located in the Uppland 
County (central-eastern Sweden). Herds were members 


of the local livestock association (Svenska Husdjur) 
and enrolled in the National Animal Disease Recording 
System (Emannelson, 1988) as well as in the Swedish 
Official Milk Recording Scheme (Olsson et ah, 2001). 
These herds participated in a specific study aiming at 
investigating the dynamics of BRSV and BCV infections 
[A. Ohlson, E. Emannelson, M. Traven, and S. Alenins 
(both of Department of Clinical Sciences, Swedish 
University of Agricultural Sciences, Uppsala, Sweden); 
unpublished data]. During this study, periodic testing 
at 4 sampling times (fall 2006, spring 2007, fall 2007, 
and spring 2008) was performed to measure the level 
of BCV- and BRSV-specific antibodies in milk. In each 
herd, information on individual milk production, SCC, 
and reproductive events (calvings, AI, and pregnancy 
checks), as well as entrance dates and exit dates and 
reasons for exit was available for all animals present at 
any point in time. 

Health Data Collection 

In each herd, at each sampling time, the level of 
BCV- and BRSV-specific antibodies was assessed from 
pooled milk samples of 5 home-bred primiparons cows. 
Testing of such samples by ELISA has recently been 
recognized as a suitable method to detect the putative 
presence of cows with antibodies to BCV or BRSV in 
dairy herds: there was good agreement between results 
of individual samples and the results of pooled samples 
(Ohlson et ah, 2009). Such a strategy has been success¬ 
fully used to determine BVDV herd status on a large 
scale (Lindberg and Alenins, 1999). Compared with 
assessment of bulk tank milk, a more immediate assess¬ 
ment of herd status was achieved by sampling young 
cows, as the antibodies may remain detectable for years 
without reinfection. The sampling was performed by 
veterinarians and technicians. Ten-milliliter test tubes 
containing 1.5 mg of the preservative agent Bronopol 
(2-bromo-2-nitropropane-1.3-diol, Steins Laboratory, 
Jonkoping, Sweden) were used. The samples were not 
diluted or centrifuged and were stored at — 20°C until 
analysis. The milk samples were analyzed for presence 
of immunoglobulin G antibodies to BCV (Alenins et 
ah, 1991) and BRSV (Elvander et ah, 1995) by com¬ 
mercially available indirect ELISA (Svanova Biotech, 
Uppsala, Sweden). The same batch was used for all 
analyses. The optical density (OD) at 450 nm was cor¬ 
rected by subtraction of the negative control antigen 
OD. To adjust for day-to-day variations, the percentage 
positivity (PP) was calculated as follows: (corrected 
OD/positive control corrected OD) x 100. Cut-off was 
set to PP equal to 20, corresponding to a corrected 
OD of 0.20, which is the cut-off recommended by the 
manufacturer for individual milk samples. Using serum, 
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Table 1 . Distribution of herd-periods 1 according to their infection status with respect to bovine respiratory 
syncytial virus (BRSV) and bovine coronavirus (BCV) 2 


BRSV infection status 1 



BCV infection status 1 


i—i 

PSI 

FREE 

RI-Y 

RI-N 

PSI-Y 

PSI-N 

h—i 

RI-Y 

1 

0 

0 

3 

1 


RI-N 

4 

1 

1 

5 

1 

PSI 

PSI-Y 

2 

0 

0 

0 

1 


PSI-N 

9 

3 

8 

26 

17 

FREE 


1 

3 

2 

3 

9 


1 See main text for definition. 

RI = presumably recently infected; RI-Y = presumably recently infected with an outbreak; RI-N = presum¬ 
ably recently infected without any outbreak; PSI = presumably past steadily infected; PSI-Y = presumably 
past steadily infected with an outbreak; PSI-N = presumably past steadily infected without any outbreak; 
FREE = presumably free of infection. 


the sensitivity and specificity of the BCV assay, rela¬ 
tive to a virus neutralization test, were estimated to be 
84.6 and 100%, respectively (91 sera from 9 Swedish 
dairy herds; Alenins et ah, 1991). Corresponding values 
for the BRSV assay were 94.0 and 100%, respectively, 
relative to an established ELISA in use at the Central 
Veterinary Laboratory (Weybridge, UK; 157 sera from 
Swedish and English herds; Elvander et ah, 1995). The 
2 assays have been recognized to work equally well in 
milk and serum [S. Alenins (Department of Clinical 
Sciences, Swedish University of Agricultural Sciences, 
Uppsala, Sweden), personal communication; Elvander 
et ah, 1995]. The same PP cut-off of 20 was applied 
to pooled milk samples, because Ohlson et al. (2009) 
found a good agreement (97 and 94%) between results 
of individual and pooled samples for BCV and BRSV 
assays, respectively. Therefore, a herd was considered 
as positive if the PP value was >20, and negative oth¬ 
erwise. 

At the end of each winter season (April 2007 and 
2008), farmers completed a questionnaire in which they 
were asked to declare the occurrence of any clinical 
signs presumably related to BCV or BRSV outbreaks, 
defined as diarrhea or respiratory signs affecting at 
least 25% of the dairy cows present within 1 wk. In case 
of occurrence, the farmers had, through not mutually 
exclusive closed questions, 1) to describe the observed 
signs (diarrhea, cough) and 2) to score their severity 
(from 0 = no to 3 = severe). 

BRSV and BCV Infection Herd Status 

Two periods-at-risk for a herd to get an infection 
were determined for each herd: fall 2006 to spring 2007 
(period A), and fall 2007 to spring 2008 (period B). 
The BRSV and BCV infection herd statuses were then 
defined on a period-at-risk basis. Therefore, in each 


herd that was sampled 4 times, 2 periods-at-risk were 
considered. 

For each herd-period-at-risk, 3 types (I, II, III) of 
herd status were defined in the same way for BCV and 
BRSV infections. These types were based on 1) the 
dynamics of BCV- and BRSV-specific antibody levels 
(type I); 2) the possible occurrence of a presumably 
BCV- or BRSV-related clinical outbreak (type II); and 
3) the combination of both pieces of information (type 
III). 

The type I herd status had 4 classes for BCV and 
BRSV infections separately. Herds that tested negative 
in fall of year n and positive in spring of year n + 
1 were presumed to be recently infected (RI); those 
that tested positive twice consecutively were presumed 
to be past steadily infected (PSI); those that tested 
positive and then tested negative were presumed to be 
past infected-recently recovered (PIRR); those that 
tested negative twice consecutively were presumed to 
be noninfected (NI), at least for the last 2 yr. Because 
our objective was to assess the production effects as¬ 
sociated with a presumably steady or recent infection, 
PIRR and NI herd-periods were considered together in 
one group taken as the reference class (FREE; Table 

a 

The type II herd-status had 2 classes separately for 
BRSV and BCV infections: absence (O-N) versus pres¬ 
ence (O-Y) of a presumably related clinical outbreak 
within the herd-period-at-risk. The information about 
the occurrence of clinical signs was used to define 
presumably BCV-related outbreaks and presumably 
BRSV-related outbreaks, according to classification 
rules, relying on the type of observed signs (diarrhea 
vs. cough) and their scoring (Figure 1). When farm¬ 
ers scored the same grade for cough and diarrhea, the 
concerned outbreaks were defined to be presumably 
BCV related, because BRSV infection is mainly associ- 
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Figure 1 . Classification of presumably bovine respiratory syncytial 
virus (BRSV)-related outbreaks (dark gray) and bovine coronavirus 
(BCV)-related outbreaks (light gray) according to the type of ob¬ 
served signs (diarrhea and cough) and their severity scoring (0 to 3). 
One herd (*) being antibody negative to BRSV and positive to BCV in 
the follow-up milk sample was considered as having a BCV outbreak. 


ated with respiratory signs, whereas BCV infection can 
cause both respiratory and enteric signs. 

The type III herd-status had 5 classes (Table 1): RI 
herd without any declared outbreak within the herd- 
period-at-risk (RI-N); RI herd with a declared out¬ 
break within the herd-period-at-risk (RI-Y); PSI herd 
without any declared outbreak within the herd-period- 
at-risk (PSI-N); PSI herd with a declared outbreak 
within the herd-period-at-risk (PSI-Y); and FREE 
herds. As expected, no FREE herds (defined based on 
the dynamics of BCV- and BRSV-specific antibody 
levels) experienced a declared outbreak. 

Modeling 

General Principle . The aim of this study was to as¬ 
sess effects associated with presumably recent or steady 
BCV and BRSV infections. Because our assumptions 
were that milk yield, SCC, reproductive failure after 
first AI, and young stock mortality may be affected by 
both infections, all models were run with both BRSV 
and BCV infection herd statuses introduced as covari¬ 
ates. Then, for BRSV infection, the provided estimates 
were adjusted for the possible effect of BCV infection 
and vice versa. For each outcome variable, 3 different 
analyses were performed: model I, model II, model III 


including type I, type II, and type III statuses defined 
above as covariates, respectively. 

The infection statuses were defined at the herd level, 
whereas the performance of cows (milk yield, SCC, 
reproduction) and of young stock (mortality) was 
considered at the individual level. Because the time of 
putative infection within a given period-at-risk was not 
known, we considered that all individual performances 
occurring within the period-at-risk were potentially af¬ 
fected. Therefore, all performances occurring between 
the starting and ending dates of the period-at-risk were 
taken into account in the analyses. Starting and ending 
dates were defined to be October 1, year n, and April 
30, year n + 1, because the fall milk samples were col¬ 
lected in September and the spring samples in May. 

The corresponding herd statuses were then assigned 
to each individual performance of each animal, depend¬ 
ing on the herd to which it belonged. Because the in¬ 
dividual performances were clustered within cow (test- 
day milk yield, test-day SCC), and the animals (cows, 
young stock) were clustered within herd, hierarchical 
mixed models considering herd and cow as random ef¬ 
fects were performed. 

Association with Milk Yield and SCC. The 

outcome variables were the individual monthly test- 
day milk yield (TDMY) and SCC (TDSCC) of cows. 
Only test-days occurring between October 1, 2006, and 
April 30, 2007, or October 1, 2007, and April 30, 2008 
were considered for analyses. The associations between 
herd-status and TDMY or TDSCC (after logarithmic 
transformation) were assessed using mixed linear mod¬ 
els, after adjustment for herd, cow, DIM, parity, breed, 
and season. The models (MIXED procedure, SAS Insti¬ 
tute Inc., Cary, NC) were as follows: 


Y = m + BRSV, + BCV, + DIM fc + LN, + BR m 


+ SEAS n + H 0 + COW„ + e 


ijklmnopi 


where Y is the TDMY or the log-transformed TD¬ 
SCC; m is the overall mean; BRSV^ is the fixed effect 
of BRSV-related herd-status i (3 classes for type I, 2 
classes for type II, 5 classes for type III); BCV, is the 
fixed effect of BCV-related herd-status j (3 classes for 
type I, 2 classes for type II, 5 classes for type III); 
DIM fc is the fixed effect of class k of DIM (13 classes); 
LN i is the fixed effect of class l of lactation number 
(5 classes: 1, 2, 3, 4, and >5); BR m is the fixed effect 
of class m of breed (3 classes: Swedish Red, Swedish 
Holstein, others); SEAS n is the fixed effect of class n 
of season (2 classes: October to December, January to 
April); H 0 is the random effect of herd o (65 herds); 
COW,, is the random effect of cow p within herd o; and 
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e ijkimno P is the residual. For each class variable, dummy 
coding was done setting the reference level to zero: each 
estimate thus measures the effect of each nonreference 
level compared with the reference level. The cow taken 
as reference was in a FREE herd when considering 
type I and III statuses, and in an O-N herd for type II 
status. To account for correlations between test days 
within cow, a compound symmetry variance matrix was 
used. A first-order autoregressive variance matrix was 
initially attempted for TDMY because other research 
(Grohn et ah, 1999) has indicated that it may be the 
most appropriate, but the models did not converge. 

Association with Reproductive Performance. 
The reproductive performance was assessed by the out¬ 
come “reproductive failure versus success.” A success 
was defined as a new calving occurring after the first 
AI after calving; that is, which was not followed by 
any other subsequent AI within the lactation and with 
the calving occurring within a delay compatible with a 
gestation length (<300 d). Cows with no information 
on subsequent calving but a successful pregnancy check 
after first AI and no subsequent recorded AI were con¬ 
sidered as experiencing a success. Other situations were 
considered as failure. However, cows with no subsequent 
events (AI, pregnancy checks, culling, calving) recorded 
after first AI or culled within 200 d after first AI were 
considered as noninformative, because their outcomes 
could not be reliably defined; they were then excluded 
from the analyses. Preliminary investigations showed 
that the proportions of noninformative animals did not 
differ statistically (P > 0.05) according to the BCV 
and BRSV infection status of the herd to which they 
belonged. Only first AI occurring between October 1, 
2006, and April 30, 2007, or between October 1, 2007, 
and April 30, 2008 were considered for analyses. The 
association between herd status and “calving after first 
service” was assessed using a mixed logistic regression 
model. The model (GLIMMIX macro, SAS Institute 
Inc.) was as follows: 

Y = m + BRSV, + BCV, + C-AI* + LN, + BR m 
+ SEAS n + MY 0 + EBp + H g + COW r + ep/ mnopgr , 

where Y is the dependent variable (failure vs. suc¬ 
cess); m is the overall mean; BRSV^ is the fixed effect 
of BRSV-related herd status i (3 classes for type I, 2 
classes for type II, 5 classes for type III); BCV, is the 
fixed effect of BCV-related herd status j (3 classes for 
type I, 2 classes for type II, 5 classes for type III); C-AI*. 
is the fixed effect of class k of calving-to-AI interval (7 
classes); LN i is the fixed effect of class l of lactation 
number (5 classes: 1, 2, 3, 4, and >5); BR m is the fixed 
effect of class m of breed (3 classes: Swedish Red, Swed¬ 


ish Holstein, others); SEAS n is the fixed effect of class 
n of season (2 classes: October to December, January 
to April); MY 0 is the fixed effect of class o of milk 
yield, expressed as the maximum of the first 3 TDMY 
(3 classes); EB,, is the fixed effect of class p of energy 
balance in early lactation, estimated as the minimum of 
the milk protein content Tat content ratio out of the first 
3 test-day records (3 classes); H, is the random effect 
of herd q (65 herds); COW r is the random effect of cow 
r within herd q\ and e^ mn(W is the residual. The cow 
taken as a reference had a success; that is, a calving or 
positive pregnancy check after first AI, in a FREE herd 
when considering type I and III statuses, and in an O-N 
herd for type II status. For each covariate, the logistic 
regression models provide adjusted odds ratios (OR). 
Because OR are overestimating relative risks (RR), 
especially when the incidence of the event under study 
is high, OR were converted into RR using the method 
developed by Beaudeau and Fonrichon (1998). 

Association with Young Stock Mortality. The 
risk of mortality in young stock was assessed using 
survival analysis to account for both the length of 
exposure to a risk of mortality that varied between 
animals, and censorship. Young stock fattened and sent 
to slaughter, sold to another herd, or temporarily out 
of the herd were considered as censored at their date 
of exit from the herd. As male calves are usually sold 
early in life in Swedish dairy herds, the young stock 
considered here were mostly female calves and heifers. 
Calves and heifers still alive at the end of the period- 
at-risk were considered as censored at its ending date 
(April 30, 2007 for period A; April 30, 2008 for period 
B). Animals born before the starting date of the period- 
at-risk (October 1, 2006 for period A; October 1, 2007 
for period B) were considered at risk of mortality in 
association with BRSV and BCV infections from this 
starting date and therefore introduced in the survival 
model as left-truncated individuals. Only calves aged 
3 to 500 d old at the ending day of the period-at-risk 
were retained for analyses, because the risk of dying of 
pneumonia or diarrhea was found to be extremely low 
afterward (Svensson et ah, 2006). 

The association between herd-status and mortality 
was assessed using a Cox proportional-hazard model 
(PHREG procedure, SAS Institute Inc.). This model 
was based on the concept of hazard function X(t), 
where \(t) is the probability for an animal of dying 
at t, given that it is still alive before t. This hazard 
characterizes the relative risk of death. The hazard rate 
\(£) = X[t, zft)] is the product of a baseline function 
X 0 (t) and of exp (<(3), representing the vector of covari¬ 
ates that influences the incidence rate of death over 
time. The effects of these covariates can be estimated 
independently from the baseline hazard function with a 
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procedure of semiparametric estimation, which involves 
the maximization of a partial likelihood. The baseline 
hazard function can then be unspecified. 

The model considered here was as follows: 

\(t) = \ Q (t) exp (BRSV, + BCVj + LN fc 

+ SEAS/) + e vjkh 

where \(t) is the hazard function at time t; X 0 (t) is the 
unspecified hazard function; BRSV, is the fixed effect 
of BRSV-related herd-status i (3 classes for type I, 2 
classes for type II, 5 classes for type III); BCV 7 is the 
fixed effect of BCV-related herd-status j (3 classes for 
type I, 2 classes for type II, 5 classes for type III); LN*. 
is the fixed effect of class k of the parity of the mother 
(2 classes: primiparous, multiparous); SEAS/is the fixed 
effect of class l of season of birth (3 classes: May to Au¬ 
gust, September to November, December to April); and 
e t j k i is the residual. The observations analyzed were the 
length of life, defined for each animal as the number of 
days between d 3 after birth and the mortality day, or if 
still alive, the last day of the period-at-risk. The assump¬ 
tion of proportional hazards was checked by including 
time-dependent covariates, defined as interactions of 
the 4 putative predictors described above and the loga¬ 
rithm of the survival time, in the Cox models and test¬ 
ing for significance. The assumption of proportionality 
was not rejected (P > 0.05) in any model (results not 
shown). To account for clustering (i.e., animals are not 
independent within herd), the data were fitted by us¬ 
ing the robust sandwich estimate (Wei et ah, 1989) for 
statistical inference (computation of robust estimates 


of standard errors). An adjusted relative hazard ratio 
(HR) was estimated for each covariate by taking the 
exponent of the estimates of effects. This HR measures 
the instantaneous hazard for a calf to die when exposed 
(e.g., in RI or PSI herds) versus being unexposed (e.g., 
in FREE herds). 

RESULTS 

A total of 101 herd-periods (55 from period A; 46 from 
period B) from 65 herds were considered for analysis. 
The different number of observations between periods 
was because of missing milk samples at some sampling 
times. The distribution of herd-periods in each BCV 
and BRSV herd status is shown in Table 1. The relative 
proportions of periods-at-risk in status RI, PSI, and 
FREE were 17, 65, and 18% and 24, 47, and 29% for 
BRSV infection and BCV infection, respectively. The 
proportion of periods-at-risk with an outbreak classi¬ 
fied as BRSV related was 8%, whereas this proportion 
was 34% for BCV related outbreaks. 

A total of 36,184 test-days, 2,716 cows with a first AI 
within the periods-at-risk, and 4,104 calves and heifers 
were retained for analyses of TDMY and TDSCC, risk 
of failure, and risk of mortality, respectively. 

Effects Associated with BRSV Infection Herd Status 

The relationships between BRSV infection herd sta¬ 
tuses and milk yield, SCC, calving after first AI and 
young stock mortality are displayed in Tables 2, 3, 4, 
and 5, respectively. 


Table 2. Adjusted effects on test-day milk yield (TDMY; kg/d) according to herd status with respect to 
bovine respiratory syncytial virus (BRSV) and bovine coronavirus (BCV) infection 




BRSV 


BCV 

Herd status 1 

TDMY 

95% Cl 

TDMY 

95% Cl 

Type I status (model I) 2 

RI 

-0.57 

-1.05; -0.08 

0.01 

-0.41; 0.42 

PSI 

0.16 

-0.36; 0.68 

-0.13 

-0.57; 0.31 

FREE 

Reference 

— 

Reference 

— 

Type II status (model II) 2 

O-Y 

-0.91 

-1.34; -0.49 

0.06 

-0.20; 0.32 

O-N 

Reference 

— 

Reference 

— 

Type III status (model III) 2 

RI-Y 

-1.42 

-2.11; -0.74 

0.32 

-0.13; 0.78 

RI-N 

-0.58 

-1.11; -0.04 

-0.95 

-1.62; -0.28 

PSI-Y 

-0.99 

-2.19; 0.21 

0.16 

-0.45; 0.76 

PSI-N 

-0.19 

-0.74; 0.36 

0.05 

-0.41; 0.50 

FREE 

Reference 

— 

Reference 

— 


RI = presumably recently infected; PSI = presumably past steadily infected; FREE = presumably free of 
infection; O-Y = with an outbreak; O-N = without any outbreak; RI-Y = presumably recently infected with 
an outbreak; RI-N = presumably recently infected without any outbreak; PSI-Y = presumably past steadily 
infected with an outbreak; PSI-N = presumably past steadily infected without any outbreak. 

See main text for definition. 
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The first, second (median), and third qnartiles of ob¬ 
served TDMY of cows in FREE herds were 2E9, 27.8, 
and 34.0 kg/d, respectively. Cows in herds classified as 
RI had a significantly (P < 0.05) lower milk yield than 
cows in FREE herds (reduction of 0.57 kg/d; Table 2). 
The reduction in milk yield in cows from O-Y herds 
(with a BRSV-related outbreak) was 0.91 kg/d (P < 
0.05). Cows in RI herds with an outbreak (RI-Y herds) 
showed a greater reduction in milk yield (—0.84 kg/d; 
P < 0.05) compared with cows in RI herds without any 
declared BRSV-related outbreak (RI-N herds). 

The geometric mean (95% Cl) of observed TDSCC of 
cows in FREE herds was 82,000 (7,000-1,020,000 cells/ 
mL). The SCC of cows in herds classified as RI and PSI 
were not significantly (P > 0.05) different from those 
in cows from FREE herds (—3,000 and 0 cells/mL, re¬ 
spectively; Table 3). However, there was a significant 
(P < 0.05) increase in SCC, corresponding to 12,000 
cells/mL, in cows from O-Y herds compared with cows 
in O-N herds. In addition, cows in RI-Y herds had sig¬ 
nificantly increased SCC (18,000 cells/mL; P < 0.05) 
compared with cows in RI-N herds. 

The proportion of reproductive failure after first AI 
of the 410 cows located in FREE herds was 60.0%. The 
BRSV-related statuses, regardless of their definition, 
were not significantly (P > 0.05) associated with calv¬ 
ing rate (Table 4). However, there was a trend toward 
a higher risk of failure in cows from herds experiencing 
an outbreak (O-Y and RI-Y) compared with O-N and 
FREE herds. 

In the population of 569 calves and heifers present in 
the FREE herds within the periods-at-risk, the mortal¬ 
ity probability at 500 d of age was 0.032 (product-limit 


estimate); that is, 3.2% of the animals at risk had died 
at that age. The BRSV-related statuses, regardless of 
their definition, were not significantly (P > 0.05) asso¬ 
ciated with the death hazard (Table 5). However, there 
was a trend toward a greater risk of death in calves 
and heifers from infected herds. Calves and heifers in a 
PSI-Y herd had a significantly higher (P < 0.05) death 
hazard (HR = 2.14) than those in a FREE herd. 

Effects Associated with BCV Infection Herd Status 

The relationships between BCV infection herd status 
and milk yield, SCC, calving after first AI, and young 
stock mortality are displayed in Tables 2, 3, 4 and 
5, respectively. The first qnartile, median, and third 
qnartile of observed TDMY of cows in FREE herds 
were 21.6, 27.4, and 34.0 kg/d, respectively. The BCV- 
related statuses, regardless of their definition, were not 
significantly (P > 0.05) associated with TDMY (Table 
2). Only cows in RI-N herds had a significant reduction 
in milk yield (—0.95 kg/d; P < 0.05) compared with 
cows in FREE herds. 

The geometric mean (95% Cl) of observed TDSCC 
of cows in FREE herds was 86,000 (7,000-1,140,000 
cells/mL). The TDSCC did not differ significantly (P 

> 0.05) between the herd statuses (Table 3). The pro¬ 
portion of failure after first AI of the 552 cows located 
in FREE herds was 56.3%. The BCV-related statuses, 
regardless of their definition, were not significantly (P 

> 0.05) associated with calving rate (Table 4). 

In the population of 833 calves and heifers present in 
the FREE herds within the periods-at-risk, the mortal¬ 
ity probability at 500 d of age was 0.048 (product-limit 


Table 3. Adjusted effects on test-day SCC (TDSCC; x 1,000 cells/mL) according to the herd-status with 
respect to bovine respiratory syncytial virus (BRSV) and bovine coronavirus (BCV) infection 




BRSV 

BCV 


Herd status 1 

TDSCC 

95% Cl 

TDSCC 

95% Cl 

Type I status (model I) 2 

RI 

-3 

-18; 4 

-1 

-13; 10 

PSI 

0 

-11; 17 

4 

-1; 28 

FREE 

Reference 

— 

Reference 

— 

Type II status (model II)“ 

O-Y 

12 

3; 34 

-3 

-19; 3 

O-N 

Reference 

— 

Reference 

— 

Type III status (model III) 2 

RI-Y 

8 

3; 52 

1 

-8; 18 

RI-N 

-10 

-34; -15 

-6 

-31; 1 

PSI-Y 

6 

-13; 65 

1 

-12; 23 

PSI-N 

0 

-13; 16 

3 

-4; 23 

FREE 

Reference 

— 

Reference 

— 


X RI = presumably recently infected; PSI = presumably past steadily infected; FREE = presumably free of 
infection; O-Y = with an outbreak; O-N = without any outbreak; RI-Y = presumably recently infected with 
an outbreak; RI-N = presumably recently infected without any outbreak; PSI-Y = presumably past steadily 
infected with an outbreak; PSI-N = presumably past steadily infected without any outbreak. 

See main text for definition. 
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Table 4. Adjusted relative risk of reproductive failure after first AI according to herd status with respect to 
bovine respiratory syncytial virus (BRSV) and bovine coronavirus (BCV) infection 



BRSV 


BCV 


Herd status 1 

Relative risk 

95% Cl 

Relative risk 

95% Cl 

Type I status (model I) 2 

RI 

1.01 

0.86; 1.21 

0.96 

0.83; 1.11 

PSI 

0.96 

0.91; 1.01 

1.04 

0.97; 1.12 

FREE 

Reference 

— 

Reference 

— 

Type II status (model II) 2 

O-Y 

1.16 

0.85; 1.58 

0.93 

0.79; 1.10 

O-N 

Reference 

— 

Reference 

— 

Type III status (model III) 2 

RI-Y 

1.07 

0.74; 1.56 

0.92 

0.76; 1.12 

RI-N 

1.02 

0.81; 1.31 

0.99 

0.68; 1.43 

PSI-Y 

0.98 

0.59; 1.67 

1.16 

0.87; 1.56 

PSI-N 

0.96 

0.91; 1.02 

1.04 

0.96; 1.14 

FREE 

Reference 

— 

Reference 

— 


X RI = presumably recently infected; PSI = presumably past steadily infected; FREE = presumably free of 
infection; O-Y = with an outbreak; O-N = without any outbreak; REY = presumably recently infected with 
an outbreak; RI-N = presumably recently infected without any outbreak; PSI-Y = presumably past steadily 
infected with an outbreak; PSI-N = presumably past steadily infected without any outbreak. 

See main text for definition. 


estimate); that is, 4.8% of the animals at risk had died 
at that age. In contrast to our assumptions, calves and 
heifers in presumably BCV-infected herds (O-Y, RI, 
and RI-Y) had a significantly lower probability of dying 
than the ones in O-N and FREE herds. 

Effects of Adjustment Variables 

The direction of estimates associated with the ad¬ 
justment covariates was comparable to that described 
in the literature. Thus, TDMY and TDSCC were lower 


in lactation number <2 than in older cows and varied 
with DIM in agreement with an expected lactation 
curve. Risk of reproductive failure was higher in cows 
in lactation number >3, with a shorter calving-to- 
AI interval, a higher milk yield, and a lower protein 
content Tat content ratio in early lactation. Test-day 
milk yield, TDSCC, and risk of reproductive failure 
were higher in winter (January to April) than in the 
period October to December. The probability of death 
was higher in calves and heifers born to primiparons 
cows. 


Table 5. Adjusted hazard ratio of calf and heifer mortality according to herd status with respect to bovine 
respiratory syncytial virus (BRSV) and bovine coronavirus (BCV) infection 



BRSV 


BCV 


Herd-status 1 

Hazard ratio 

95% Cl 

Hazard ratio 

95% Cl 

Type I status (model I) 2 

RI 

1.42 

0.73; 2.77 

0.50 

0.26; 0.96 

PSI 

1.26 

0.77; 2.04 

1.19 

0.69; 2.05 

FREE 

Reference 

— 

Reference 

— 

Type II status (model II) 2 

O-Y 

1.53 

0.85; 2.76 

0.45 

0.25; 0.80 

O-N 

Reference 

— 

Reference 

— 

Type III status (model III) 2 

RI-Y 

1.58 

0.63; 3.95 

0.50 

0.25; 1.00 

RI-N 

1.19 

0.54; 2.62 

0.39 

0.14; 1.11 

PSI-Y 

2.14 

1.29; 3.54 

0.59 

0.21; 1.66 

PSI-N 

1.18 

0.77; 1.80 

1.37 

0.79; 2.39 

FREE 

Reference 

— 

Reference 

— 


X RI = presumably recently infected; PSI = presumably past steadily infected; FREE = presumably free of 
infection; O-Y = with an outbreak; O-N = without any outbreak; RI-Y = presumably recently infected with 
an outbreak; RI-N = presumably recently infected without any outbreak; PSI-Y = presumably past steadily 
infected with an outbreak; PSI-N = presumably past steadily infected without any outbreak. 

See main text for definition. 
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DISCUSSION 

This study identified losses in performance associated 
with BRSV infection, whereas our estimates regarding 
BCV-infected herds were either not significantly dif¬ 
ferent from those in noninfected herds or contrary in 
direction to biological expectations. 

A significant reduction in milk yield was observed 
for cows in presumably BRSV RI herds, as well as in 
presumably BRSV O-Y herds compared with cows 
in presumably infection-free herds (FREE and O-N 
herds, respectively). In addition, there was a further 
significant detrimental effect in presumably RI herds 
when they experienced an outbreak (Table 3). These 
results are in accordance in direction with those re¬ 
ported by Norstrom et al. (2001) in case of a declared 
outbreak, but the magnitude of the decrease observed 
here (on average, a 0.9 kg/d reduction over a 7-mo 
period) was much higher than that in the latter study 
(0.7 kg/d over a 1-wk period). However, Norstrom et al. 
(2001) hypothesized that their result could be largely 
underestimated; their estimate was similar to the one 
reported by Van der Poel et al. (1995) in case of BRSV 
reinfections. Information about the relative proportion 
of primary infected and reinfected cows (assumed to be 
severely affected and to have a smaller loss in milk yield, 
respectively) in the affected herds was lacking in their 
study. The results of our study could also support their 
assumption, as the milk yield of cows in presumably 
past steadily infected herds (comprising either past- 
infected cows or to some extent reinfected cows) was 
not different from that in presumably FREE herds. 

The risk of reproductive failure after first service 
(Table 4) as well as risk of death in calves and heifers 
(Table 5) was found to be increased in herds with a pre¬ 
sumably BRSV-related outbreak, for the former, and 
in presumably recently or past-infected herds, for the 
latter. These findings were in agreement with biologi¬ 
cal expectations and with the literature. The mortality 
rate can reach 20% in some severe outbreaks (Stott et 
ah, 1980), and pneumonia and cough are recognized as 
major risk factors for death early in life (Waltner-Toews 
et al., 1986; Curtis et ah, 1988; van der Fels-Klerx et 
ah, 2002; Svensson et al., 2006). These estimates were, 
however, not significantly different from unity, possibly 
because of a lack of statistical power because of the 
relatively low sample sizes (2,716 cows with a first Al 
within the periods-at-risk and 4,104 calves and heifers 
at risk of mortality, respectively). Further studies are 
needed to confirm or reject the trends observed here. 
We also observed a significant increase in SCC (of more 
than 10,000 cells/mL) for cows in herds with a BRSV 
outbreak. This increase cannot only be caused by a 
concentration effect, because the increase would have 


been about 3,000 cells/mL if exclusively in relation to 
the associated reduction in milk yield (—0.91 kg/d) 
observed in the same herds. 

The present study showed that, regardless of the 
performance under study, the point estimates of effects 
associated with BRSV infection were large in magni¬ 
tude, despite the fact that they were calculated based 
on animal performance in the whole period-at-risk (7 
mo duration). This finding, which was not expected 
because the BRSV infection was previously found to 
have only a temporary effect (Norstrom et al., 2001), 
suggests that the effect of this infection may be very 
severe at the time of its occurrence. 

It should be discussed, however, whether the losses in 
productivity estimated in this study could be regarded 
as a direct consequence of the BRSV infection only, be¬ 
cause causal relationships cannot easily be established 
in observational studies such as this. One possibility is 
that the increase in SCC, the higher risk of failure after 
first service, and higher risk of mortality of young stock 
are markers of poor hygiene and management. Poorly 
managed herds might also lack biosecurity and are 
therefore more likely to get infectious diseases. How¬ 
ever, all herds under study were free from bovine viral 
diarrhea virus, bovine paratuberculosis, and infectious 
bovine rhinotracheitis. 

The same question could be raised regarding BCV 
infection. However, the present study provided either 
nonsignificant estimates (e.g., for milk yield and SCC) 
or significant estimates that were contrary in direction 
to biological hypotheses and knowledge (e.g., lower risk 
of young stock death in presumably RI and in herds 
with a presumably BCV-related clinical outbreak; 
Table 5). 

Two explanations could be put forth for the absence 
of any observed decrease in milk yield associated 
with BCV, although a decrease could be expected 
from farmers’ experiences as well as from description 
of BCV-related outbreaks in the literature. First, the 
decreased milk yield associated with winter dysentery, 
even if large in magnitude (on average 22%), is often 
described as being of short duration (<2 wk; Traven, 
2000). Therefore, if it exists, such an effect would have 
been difficult to detect using test-day records, which 
are usually performed 1 mo apart. Second, the TDMY 
used for comparisons were all those within the period- 
at-risk (i.e., from October year n to April year n + 1), 
whereas winter is known to be a predisposing factor for 
BCV outbreaks (Saif, 1990). Thus, this strategy may 
have “smoothed” the possible decrease in milk yield 
associated with an outbreak over the whole period-at- 
risk. However, no precise (and reliable) information on 
actual dates of possible outbreaks were available and 
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a more detailed analysis could therefore not be per¬ 
formed. 

This second explanation may be less likely, because 
the milk yield in presumably RI-N herds was found to 
be significantly lower than in RI-Y (Table 2). Because 
both these groups experienced a recent BCV circula¬ 
tion, this finding might suggest that herds in which the 
owner noticed a presumably BCV-related clinical out¬ 
break are better managed to some extent with higher 
performance. This hypothesis is also supported by the 
fact that only 3 of the 28 herds with a presumably 
BCV-related clinical outbreak also had a presumably 
BRSV-related clinical outbreak (for which a decrease in 
milk yield was observed; Table 1). In addition, it may 
be questioned if all outbreaks classified as BCV related 
were indeed related to BCV infection, because signs 
of diarrhea (the main criterion to define an outbreak) 
may well be related to other factors; for instance, re¬ 
lated to feeding regimens, especially in high-performing 
herds where cows are more sensitive. This reason for 
misclassification was deemed not differential between 
herds exhibiting or not exhibiting specific antibodies 
to BCV, because the farmers were not aware of the 
status of their herd when they declared the observed 
clinical signs. In that context, the presumably BCV- 
related outbreaks might have been overestimated in our 
study. The same may not necessarily apply for cough, 
for which the risk of misclassification toward BRSV 
infection was deemed lower. 

The assumption of better management and perfor¬ 
mance in herds with a presumably BCV-related out¬ 
break could also indirectly explain, at least partially, 
the trend toward a lower risk of death in young stock in 
such herds (Table 5). This might also be a reason why, 
contrary to our expectations, BCV outbreaks were not 
found to be associated with increased SCC; the higher 
level of management might prevent the higher risk for 
a cow with diarrhea to be exposed to pathogens likely 
to infect the udder. 

According to this study, BRSV infection had a 
large significant effect on milk yield and SCC and to 
some extent on reproduction performance of cows and 
survival of young stock. However, for the latter, the 
relationships found here should be confirmed in a larger 
study. The quantitative estimates of production effects 
associated with this infection will be useful information 
to assess its economic impact, as well as the profitabil¬ 
ity of control schemes that could be implemented on 
farms. In contrast, BCV infection, identified here based 
on either the dynamics of BCV-specific antibody levels 
over time or the occurrence of a presumably related 
outbreak, was not significantly associated with any 
production losses in cows from the infected herds com¬ 


pared with those in the FREE herds. These findings 
were unexpected because severe decrease in milk yield 
in cows with winter dysentery, as well as increased calf 
mortality, were identified in previous descriptions of 
outbreaks. Because the consequences on performances 
associated with BCV-related outbreaks are probably of 
short duration in most situations, case-control investi¬ 
gations relying on very reliable identification and dating 
of the BCV-related outbreak and of the corresponding 
period-at-risk, and targeting daily performance records 
would be a profitable area of research. 
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